Electronic Griffiths phase of the d = 2 Mott transition 
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We investigate the effects of disorder within the T = Brinkman-Rice (BR) scenario for the Mott 
metal-insulator transition (MIT) in two dimensions (2d). For sufficiently weak disorder the transition 
retains the Mott character, as signaled by the vanishing of the local quasiparticles (QP) weights Zi 
and strong disorder screening at criticality. In contrast to the behavior in high dimensions, here 
the local spatial fluctuations of QP parameters are strongly enhanced in the critical regime, with a 
distribution function P(Z) ~ Z"^ 1 and a — > at the transition. This behavior indicates a robust 
emergence of an electronic Griffiths phase preceding the MIT, in a fashion surprisingly reminiscent 
of the "Infinite Randomness Fixed Point" scenario for disordered quantum magnets. 
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The effects of disorder on quantum criticality [l| prove 
to be much more dramatic than in classical systems. 
Here, some critical points can be described by an "infinite 
randomness fixed point" (IRFP) and the associated 
quantum Griffiths phase. Such exotic behavior is well 
established in insulating quantum magnets with discrete 
internal symmetry of the order parameter [3], but may 
or may not survive in other models or in the presence of 
dissipation due to conduction electrons. 

More general insight in the robustness of the IRFP 
scenario rests on a recently proposed symmetry classifi- 
cation [l[ , based on the lower critical dimension of droplet 
excitations. These ideas have found support in very re- 
cent work sparking considerable renewed interest [H|- 
Much of this progress, however, relies on the ability to 
identify an appropriate order parameter, describing the 
corresponding symmetry breaking transitions. 

The metal-insulator transition (MIT) represents an- 
other important class of quantum criticality, one that of- 
ten cannot be reduced to breaking any static symmetry. 
Conventional theories of the MIT in disordered systems 
[fj], based on the diffusion mode picture, strongly resem- 
ble standard critical phenomena and thus do not easily 
allow [3] for rare event physics or IRFP behavior. There 
currently exists, however, a large body of experimental 
work [7j , documenting disorder-induced non- Fermi liquid 
behavior due to rare disorder configurations, even in sys- 
tems far from any spin or charge ordering. 

Theoretically, such "electronic Griffiths phases" (EGP) 
0, 0, E3] have recently been proposed for correlated elec- 
tronic systems with disorder, based on generalized dy- 
namical mean- field theory (DMFT) approaches @, ED]. 
All these works were performed on the Bethe lattice and 
identified EGPs only in the vicinity of disorder-driven 
MITs, in particular, only for strong enough disorder, in 
contrast to quantum magnets where even weak disorder 
often results in IRFP behavior. Some key unanswered 
questions thus remain: (a) What is the effect of weak 
to moderate disorder on interaction-driven MITs such as 



the Mott transition in finite dimensions? (b) Is the criti- 
cal behavior dramatically changed as in the examples of 
IRFP or may a more conventional scenario suffice? 
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Figure 1: T = phase diagram of the disordered half-filled 
Hubbard model in d = 2, as a function of the interaction U at 
weak to moderate disorder strength W < U. An intermediate 
Electronic Griffiths Phase emerges separating the disordered 
Fermi liquid metal and the Mott insulator. The inset shows 
the typical (Z tyv ) and average (Z av ) values of the local quasi- 
particle weight Zi as a function of U. The Mott transition 
is identified by the (linear) vanishing of Z typ . Note that Z av 
is finite at U c , indicating that a fraction of the sites remains 
nearly empty or doubly occupied. 

In this Letter, we investigate the effects of weak and 
moderate disorder on the Mott MIT at half filling [H] 
in two dimensions. As the simplest description of the ef- 
fects of disorder on the Mott transition, we work within 
a Brinkman-Rice (BR) scenario [ill], where a Gutzwiller 
variational approximation is applied to a disordered two 
dimensional Hubbard model. Our results demonstrate 
that: (i) for sufficiently weak disorder the transition re- 
tains the second order Mott character, where electrons 
gradually turn into localized magnetic moments; (ii) 
disorder-induced spatial inhomogeneities give rise to an 
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intermediate EGP that displays IRFP character at criti- 
cality, even when the transition is approached by increas- 
ing the interaction at weak disorder; (iii) the renormal- 
ized disorder seen by quasiparticles is strongly screened 
only at low energies, resulting in pronounced energy- 
resolved inhomogeneity of local spectral functions. 

Model. — We focus on the paramagnetic disordered 
Hubbard model with nearest-neighbor hopping and with 
site energies Ei uniformly distributed in the interval 
[-W/2,W/2] 0,[3|. We approach the Mott transition 
by increasing the on-site Hubbard interaction U at half 
filling (chemical potential [i = U/2), on an LxL square 
lattice with periodic boundary conditions. All energies 
will be expressed in units of the clean Fermi energy (half- 
bandwidth) Ep = it, where t is the hopping amplitude. 

Within our disordered BR approach, we self- 
consistently calculate the local single-particle self- 
energies Ei(w„) 
form 



id [lj] , which assume a site-dependent 
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The renormalized site energies = Vi{ei,di) and the 
local quasiparticle (QP) weights Zj = Zi (e,, df) are vari- 
ational^ calculated through the saddle-point solution of 
the corresponding Kotliar-Ruckenstein (KR) slave boson 
functional 

F = -2T^Tr\n \-iu n l + Zv + Z^HqZ 1 / 2 

+ E [ ud * - i 1 - e * + ^ - £ *+ rf] • ( 2 ) 

i 

Here, and di are the KR slave boson amplitudes flj| . 
T is the temperature, and to„ are the Matsubara frequen- 
cies. The operators Z and v are site-diagonal matrices 
[Z]ij = ZiSij; [v]ij — ViSij, and Ho is the clean and non- 
interacting lattice Hamiltonian. 

This approach is mathematically equivalent to a gen- 
eralization of the dynamical mean field theory (DMFT) 
IH to finite dimensions, the "statistical DMFT"^ im- 
plemented using a slave boson impurity solver, which pro- 
vides an elegant and efficient computational approach, 
allowing us, for example, to calculate Z t values spanning 
eight orders of magnitude. We considered several lattice 
sizes ranging up to L = 50, and for every (U, W) pair 
we typically generated around forty realizations of disor- 
der. We carefully verified that for such large lattices, all 
our results are robust and essentially independent of the 
system size (see, e.g., the inset of Fig. [2]). 

Phase diagram and Griffiths phase. — To characterize 
the T = disordered Mott transition in d — 2, we fol- 
low the evolution of the local QP weights Zi , as the in- 
teraction U is increased at fixed disorder W. For weak 
to moderate disorder, we find behavior partly reminis- 
cent of that previously established for high dimensions 
("DMFT limit") [Hj]. The approach to the critical point 



at U = U C (W) is identified by the vanishing of the typical 
QP weight Z typ = exp{(ln Z,-)}(Io|. indicating the Mott 
transmutation of a finite (large) fraction of electrons into 
local magnetic moments. Because random site energies 
tend to push the local occupation away from half filling 
U C (W) increases with disorder (Fig. [1]). 
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Figure 2: Distribution of z; = Zi/Zo (see text) as the transi- 
tion is approached by increasing U (thin lines correspond to 
U/U C {W) = 0.6,0.8,0.9,0.92,0.94,0.97). For reference, the 
thick solid line shows the DMFT "fixed point" distribution 
which remains bounded from below. Our d — 2 results show 
that, due to rare events, a low-z tail P (z) ~ emerges for 
z < 1 as the transition is approached. In the critical region 
the distribution assumes a singular form (a < 1), indicat- 
ing the onset of an Electronic Griffiths Phase. Results are 
shown for W = 5.0 and L = 20. The inset illustrates how 
for such large lattices our results for P (Z) are essentially in- 
dependent of the system size (shown for U/U C {W) = 0.94 
and a = 0.77 ± 0.06, 0.81 ± 0.05, 0.82 ± 0.05 corresponding to 
L = 20, 30, 50, respectively). 

The role of fluctuation effects, however, is best seen by 
contrasting our d = 2 results to those found in the DMFT 



limit There each site has many neighbors, and thus 
"sees" the same (self-averaged) environment ("cavity") , so 

Z DMFT = Z DMFT depends only on the local site en _ 

ergy e,. Its minimum value corresponds to the sites clos- 
est to half-filling Z™ FT ee Z = Z DMFT {e, = 0). In 
the critical region all Z? MFT - (Z) - Uf )MFT {W) - U, 
but the scaled local QP weights z f MFT = Zf MFT jZ 
approach finite values at the transition, with z^^ FT = 1. 
The corresponding scaled distribution P (zf MFT ) ap- 
proaches a fixed-point form close to U® MFT (W) (shown 
by the thick solid line in Fig. [2|. 

In low dimensions, site-to-site cavity fluctuations give 
rise to a low-Z tail emerging below the DMFT minimum 
value Z a (Fig. [2]). To bring this out, we present our 
d = 2 results in precisely the same fashion as in the 
DMFT limit, i.e. scaling each Zi with Z . Away from 
the transition the distribution resembles the DMFT form, 
but in the critical region the low-z tail assumes a power- 
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law form 



P[z) 



(3) 



Physically, the emergence of a broad distribution of 
local QP weights Zj indicates the presence of rare dis- 
order configurations characterized by anomalously low 
local energy scales Aj » Z^Ep. Since the approach to 
the Mott insulator corresponds to Z — ► 0, such regions 
with Zj <C Z typ should be recognized as "almost local- 
ized" Mott droplets. Within our BR picture, each local 
region provides a contribution Xi ~ 7i ~ to 
the spin susceptibility or the Sommerfeld coefficient, re- 
spectively. The local regions with the smallest Zj thus 
dominate the thermodynamic response and produce non- 
Fermi liauid metallic behavior [3. Iflf whenever a < 1. 
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Figure 3: Exponent a (U) for two different values of disorder. 
To calculate a we use two methods (open and closed symbols, 
see text). Result are shown for L = 20. Within the estimated 
error bars, we find an extrapolated value consistent with a = 
at the critical point U — U C (W). 

IRFP-like behavior. — To carefully calculate the expo- 
nent a(U) as the transition is approached, we use two 
distinct methods. The first relies on the "estimator" [ijj 
a = (In [Zmax/ZiDz^z^, wnere Z max ~ Z typ ~ Z a 
(see Fig. [2J is an appropriate upper-bound on the power- 
law behavior. The second approach consists in calculat- 
ing (Z _M ) z <z for given as a function of U. This 
quantity is expected to vanish at U = Z7 M , satisfying 
a {Ufj.) = 11. Both methods give consistent results (open 
and closed symbols, respectively, in Fig. [3]), which agree 
within the estimated error bars. 

The exponent a is found to decrease smoothly as the 
transition is approached, until the distribution assumes 
a singular form (a < 1), indicating the emergence of an 
Electronic Griffiths Phase (EGP) @, Its estimated 
onset (a = 1) is generally found to strictly precede the 
MIT (dashed line in Fig. [TJ, indicating that disorder fluc- 
tuations qualitatively modify the critical behavior even 
for weak to moderate disorder. Remarkably, we find that, 
within our numerical accuracy, a — > precisely along the 



critical line U = U C (W)\ This establishes a phenomenol- 
ogy which closely parallels the behavior of magnetic Grif- 
fiths phases with IRFP behavior [J, 0] . 
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Figure 4: (a) Spatial distribution of the (normalized) local 
spin susceptibilities %i ~ %i j illustrating a typical disorder 
realization containing a "rare event" (RE) with \i ^ Xtyp'-, (t>) 
Disorder fluctuations are eliminated outside a box of size I = 
9, without appreciably affecting the RE; (c) When the box is 
further reduced (here I = 3) the RE is rapidly (exponentially) 
suppressed, establishing the non-local nature of the rare event, 
in strong support of the IRFP picture. Results shown for 
L = 30, W = 5.0 and U/U c = 0.96. 



Structure of the rare events. — To explore the nature of 
the rare events (REs) dominating the EGP, i.e. the re- 
gions with Zj <C Z typ , we examined a number of disorder 
realizations and selected those few samples containing 
the smallest Zj. A typical example is shown in Fig. 01 
where the RE is seen as a very sharp peak of the local 
spin susceptibility Xi ~ Z" 1 (note the logarithmic scale). 
The corresponding RE site is then placed in the middle 
of a box of side I (dashed line in Fig. [4]). To examine the 
spatial correlations, we preserve the same disorder real- 
ization within this box, while the outside is replaced by 
an appropriate DMFT effective medium. We then recal- 
culate the QP parameters Zj as the box size is reduced 
from I — L (original model), down to I — 1 (DMFT limit 
where all spatial correlations are suppressed). We find 
that the RE is essentially unmodified until the box size 
reaches I ~ Ire (~ 9 for the example in Fig. 4) , and then 
is rapidly (exponentially) suppressed for I < Ire ■ We also 
find that the variance of the disorder strength within the 
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box of size Ire is appreciably weaker than on the average, 
establishing that the REs dominating the Griffiths phase 
stem out from rare disorder configurations, precisely as 
expected within the IRFP scenario. 

Critical behavior of the spatial inhomogeneity. — As 
a — » in the critical region, the P{Zi) distribution 
becomes "infinitely broad", since a -1 measures 0] the 
variance of lnZ^. The thermodynamic response be- 
comes increasingly inhomogeneous as the transition is 
approached; such behavior is typically seen in NMR ex- 
periments on materials displaying disorder-driven NFL 
behavior 0]. 




Figure 5: Spatial distribution of the LDOS normalized by 
its clean value, for one disorder realization, shown: (a) away 
from the Fermi energy (ui = 0.10); (b) at the Fermi energy 
(lo — 0). Due to static disorder screening, the distribution be- 
comes homogeneous close to the Fermi level, but displays pro- 
nounced spatial structures at higher energies. Results shown 
for L = 50, U/U c = 0.96 and W = 0.75. 



giving rise to a very small (but finite) renormalized dis- 
order strength W e ff = \J (vf) <C 1 at criticality (e.g. for 
W = 5, W e ft w 0.05). Near the Fermi energy {u> = 0), we 
predict the LDOS spectra to appear increasingly homoge- 
neous in the critical region. At higher energies, however, 
the very broad distribution of local QP weights (essen- 
tially local QP bandwidths) creates a very strong effec- 
tive disorder seen by the quasiparticles, and we expect 
the system to appear more and more inhomogeneous as 
criticality is approached. This result is illustrated by ex- 
plicit computation of the DOS profile (Fig. [5]), which 
is surprisingly reminiscent of recent spectroscopic images 
on doped cuprates [la ]. Our theory, which does not in- 
clude any physics associated with superconducting pair- 
ing, strongly suggests that such energy-resolved inhomo- 
geneity is a robust and general feature of disordered Mott 
systems. 

Conclusions. -We presented the first detailed model 
calculation investigating the effects of moderate disorder 
on the Mott metal insulator in two dimensions. Our find- 
ings indicate that rare disorder fluctuations may domi- 
nate quantum criticality even in absence of magnetic or- 
dering - an idea that begs experimental tests on a broad 
class of materials. The Brinkman-Rice scenario we con- 
sidered, which focuses on local (Kondo-like) effects of 
strong correlation (while neglecting inter-site magnetic 
correlations), may be relevant only for systems with suf- 
ficiently strong magnetic frustration, such as 3 He mono- 
layers adsorbed on graphite [l9| . Such variational ap- 
proach should be generalized for systems, such as copper 
oxides, where the inter-site super-exchange is strong, but 
this fascinating research direction remains a challenge for 
future work. 
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But what to expect from STM experiments directly 
measuring the local electronic spectra? Within our BR 
approach, the local density of states (LDOS) in question 
Pi{id) = — lmGii(uj — i0 + ) depends not only on the local 
QP weights Zi , but also on the renormalized site energies 
Vi through 



(Z-^-v-Hq) 



(4) 



The quasiparticles thus "see" a frequency-dependent ef- 
fective disorder potential 



e1 rs {u>)=Vi-u/Zi 



(5) 



Our explicit calculations find that the renormalized site 
energies u, become strongly screened near the transition, 
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